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PERTURBATION OF TRANSPORTATION POLYTOPES 
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Abstract. We describe a perturbation method that can be used to reduce 
the problem of finding the multivariate generating function (MGF) of a non- 
$H ' simple polytope to computing the MGF of simple polytopes. We then con- 

^ ' struct a perturbation that works for any transportation polytope. We apply 

this perturbation to the family of central transportation polytopes of order 
kn X n, and obtain formulas for the MGFs of the feasible cone of each vertex 
1/^ ' of the polytope and the MGF of the polytope. The formulas we obtain are 

enumerated by combinatorial objects. A special case of the formulas recov- 
_ ers the results on Birkhoff polytopes given by the author and De Loera and 

C J ■ Yoshida. We also recover the formula for the number of maximum vertices of 

r \ 1 transportation polytopes of order kn X n. 

^ . 

1. Introduction 

Counting number of lattice points inside polytopes is a problem of much interest. 
In Section 5 of [2] , the authors describe a method of using residue calculations to 
^ . obtain the number of lattice points of a rational polytope P from its multivariate 

'^ ' generating function or MGF f{P,z) (see Definition 2.1). In [10], the authors ex- 

tend the residue calculation described in [2] to dilations of integral polytopes P, 

f^ ' and describe how one can obtain the Ehrhart polynomial and the volume of P 

from its MGF f{P,z) (see Lemma 2.15). Therefore, the multivariate generating 
function is an interesting and useful object to study. The most common way of 

^^ . computing the MGF of a polytope (or a polyhedron) P is to compute the MGF of 

the tangent/feasible cone of each vertex of P and then sum them up (see Theorem 
2.10 and Corollary 2.11). However, computations involving non-simple cones are 
often complicated. When a cone is non-simple, a common method used to compute 
its MGF involves triangulating the cone into simple cones [2, 1]. Perturbation is 
another method often used to deal with non-simple cones in computational math- 
ematics [5, 7, 9], but it has not been widely used for computing MGFs. (However, 
see [4, 12] for some work in which perturbation is used for computing MGFs.) The 
following theorem describes a framework for globally perturbing polyhedra, with 
which wc can compute the MGFs of non-simple polyhedra. Although the theorem 
is equivalent to techniques known to the experts in the field, the particular pre- 
sentation is well suited for our purposes, and we are not aware of it being in the 
literature in this form. 

Theorem 1.1. Let A be an N x D matrix and b e K.'^. Suppose P is a non-empty 
integral polyhedron in M.^ defined by Ax < b and h{t) is a continuous function 
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on some interval containing such that b(0) = b. Suppose for each t ^ in 
the interval, h{t) defines a non-empty polyhedron P{t) = {x | Ax. < h{t)} with 
exactly £ vertices: Wt,i, ■ ■ ■ ,Wt/, and the feasible cone of P{t) at Wtj does not 
depend on t, that is, for each j '■ 1 < j < i-, there exists a fixed cone Kj such that 
iconc{P{t),wt,j) = A'j for all t ^ 0. Then 



i=i 

Note that the convergence of Wtj is proved in Theorem 3.1, which also provides 
a local version for the perturbation method. According to Theorem 1.1, if we are 
given a non-simple integral polytope P and are able to find a way of perturbing P 
into simple polytopes satisfying conditions of Theorem 1.1, we can use the MGF of 
the feasible cones of the perturbed polytope P{t) to calculate the MGF of P. In [10], 
the authors find a combinatorial expression for the MGF of the Birkhoff polytope 
Bn, from which they obtain the first combinatorial formulas for the volume and 
Ehrhart polynomial of -B„ by residue calculations. In the present paper, we apply 
our perturbation method to generalize the results in [10]. Before stating our results, 
we recall some definitions. 

The Birkhoff polytope belongs to the family of transportation polytopes: Given 
r = (ri, . . . , Tm) and c = (ci, . . . , c„) two vectors of positive entries whose coor- 
dinates sum to a fixed number, the transportation polytope determined by r and 
c, denoted by T(r,c), is the set of all m x n nonnegative matrices in which row i 
has sum r^ and column j has sum Cj. We call 7"(r, c) a transportation polytope 
of order m x n. The problem of finding MGFs of transportation polytopes is rela- 
tively easy when a transportation polytope is non-degenerate (see Section 4 for the 
definition), in which case one can apply Corollary 4.8 to find its MGF. However, if 
a transportation polytope is degenerate, it is usually non-simple, in which case one 
often has to triangulate the non-simple feasible cones. We describe a perturbation 
(in Lemma 5.1) for all the transportation polytopes so that the problem is reduced 
to finding the MGF of non-degenerate cases. In addition to this, our perturbation 
provides families of transportation polytopes that obtain the maximum possible 
number of vertices among all the transportation polytopes of the same order (see 
Lemma 5.8). 

Addressing a question of Bernd Sturmfels. we consider a special family of trans- 
portation polytopes that contains the Birkhoff polytope: the family of central trans- 
portation polytopes. A classical central transportation polytope of order to x n is 
a transportation polytope whose column sums are all m and row sums are all n. 
However, strictly speaking Birkhoff polytopes do not belong to this family because 
the column sums and row sums are all 1 for i?„. Therefore, we slightly generalize 
the definition of the classical central transportation polytopes to include Birkhoff 
polytopes. A transportation polytope T(r, c) is central of order 7ti x n if all the 
column sums are the same and all the row sums are the same. The family of central 
transportation polytopes is an interesting subset of transportation polytopes. For 
example, when m and n are coprime, the central transportation polytope of order 
m x n achieves the maximum possible number of vertices among all the transporta- 
tion polytopes of order mxn [3]. In [10], the combinatorial data used to enumerate 
the MGF of B„ is the family of rooted trees. Sturmfels asked whether we can give 
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a nice description in terms of trees for the MGF of any central transportation poly- 
tope. We answer his question for central transportation polytopes of order kn x n 
by using the perturbation method we develop. 

Denote by STk,n the set of spanning trees of the complete bipartite graph Kkn.n 
with right degree sequence (fc + 1, fc + 1, . . . , fc + 1, fc). (See Definition 5.4 for the 
definition of right degree sequence.) 

Theorem 1.2. Assume that T(r,c) is a central transportation polytope of order 
m X 71, where m = kn, and r ~ (a, . . . , a) and c = (&,..., 6) are two integer vectors. 
Suppose M is a vertex o/T(r,c). Then the MGF of the feasible cone o/7'(r, c) 
at M is 

(1.1) /(fcOne(r(r,c),Af),z)= ^ [] l _ ^cycIe(T.e) • 

TePcrtAux(A/) e<^E(T) 

Thus, the MGF of T{r, c) is 

(1.2) /(r(r,c),z)= Y: -"^ n 1 _ ,cyc.e(T,e) ■ 

T£STk,„ e<^E(T) 

In both equations, cycle(T, e) denotes the m x n {0, —1, l)-niatrix associated to the 
unique cycle in the graph TUe (see Definition 4-5 for details). In the first equation, 

Pert Aux(A/) = {T G STk,n \ the auxiliary graph of M is a subgraph of T , } 

and in the second equation, Mt is the unique vertex of T(r, c) whose auxiliary 
graph is a subgraph ofT. See Definition 4-2 for the definition of auxiliary graph. 

Note that the Birkhoff polytope i?„ is a special case of this family of central 
transportation polytopes when fc = 1. Our result recovers the formulas for the 
MGFof B„ given in [10]. 

We remark that in [11], the authors also provided a framework for global per- 
turbation method using the idea of lattice ideals (see [11, Corollary 7.5]). The 
fundamentals of these two theories are similar; we just use different approaches. 
However, in our paper, we also provide a local version for the perturbation method 
in Theorem 3.1. 

This paper is organized as follows. In Section 2, we give background on results 
related to multivariate generating functions. In Section 3, we give the framework 
for our perturbation method. In Section 4, we review results on properties of trans- 
portation polytopes. In Section 5, we give a perturbation that works for every 
transportation polytope. In Section 6, we apply the perturbation method to cen- 
tral transportation polytopes of order kn x n and give combinatorial formulas for 
the MGFs, Ehrhart polynomials, and volumes of these polytopes. When we spe- 
cialize these formulas to Birkhoff polytopes, we recover the results in [10]. We also 
recover the formula for the maximum possible number of vertices of transportation 
polytopes of order kn x n. 

2. Background 

A polyhedron is the set of points defined by a system of linear inequalities 

n,(x) < 6,;, \fl<i<N, 

where n^ is a D- vector and ni(-) is the linear function that maps x G M.^ to the dot 
product of rii and x. Since we can consider ni(-) as a point in the dual space (R")* 
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of K" , we will use the notation n; (or any bold letter) to denote both the D-vector 
and the linear function. For simplicity, we let A be the N x D matrix whose row 
vectors arc n^'s and b = (5i, . . . , bpf)'^ , so the system of linear inequalities can be 
represented as 

Ax < b. 

A polytope is a bounded polyhedron. We assume familiarity with basic definitions 
of polyhedra and polytopcs as presented in [8, 14]. For any polyhedron P, wc use 
Vert(P) to denote the vertex set of P. An integral polyhedron is a polyhedron 
whose vertices arc all lattice points, i.e., points with integer coordinates. A rational 
polyhedron is a polyhedron whose vertices arc all rational points. 

2.1. Ehrhart polynomials and multivariate generating functions. For any 

polytope P C K^, and a nonncgativc integer t, we define 

i(p,t) = #(tPnz-°) 

to be the number of lattice points inside tP = {tx \ x e P}, the tth dilation of 
P. It is well-known that given a d-dimensional integral polytope P, the function 
i(P, t) is a polynomial in t of degree d with leading coefficient being the normalized 
volume of P. Since this was first discovered by Ehrhart [6], wc often refer to i(P, t) 
as the Ehrhart polynomial of P. Because the leading coefficient of i{P, t) gives the 
volume of P, obtaining the Ehrhart polynomials of polytopcs is one way people use 
to compute volumes of polytopcs. One can find the Ehrhart polynomial i{P,t) of 
P using the multivariate generating function. 

Definition 2.1. Let P C M^ be a polyhedron. The multivariate generating func- 
tion (or MGF) of P is: 

/(P,z)= J2 ^"' 

where z" = n^i2r-. 

Note that when P is a polytope, we obtain i(P, t) by plugging Zi = 1 for all i in 
the MGF f{tP, z) of the tth dilation tP of P. 

2.2. Cones of polyhedra. One benefit of computing the MGF of a polyhedron 
is that the problem can be reduced to computing the MGF of the tangent/feasible 
cones of the given polyhedron by applying Brion's theorem. Let's first review some 
related definitions and results. 

Definition 2.2. Suppose P is a polyhedron defined by Ax < b, where the row 
vectors of A are ni, . . . , uat. Let / be a face of P. If ni(x) = hi for any x e /, 
we say the inequality ni(x) < bi supports the face / or ni(x) = 6^ is a supporting 
hyperplane for /. 
We denote by 

N(v4,b,/) = {n, |n,(x) = 6„Vxe/} 

the set of the normal vectors n^ of all the inequalities among Ax < b that support 
the face /. 

Definition 2.3. Let K C M^ be a cone. The polar cone of K is the cone 

K° = {ye{R''r |y(x)<0,VxeX}. 
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Definition 2.4. Suppose P is a polyhedron and w G P. The tangent cone of P at 

V is 

tcone(P, v) = {v + u \ v + 5u G P ioi all sufficiently small S > 0}. 
The feasible cone of P at v is 

fconc(P, w) = {u I u + (5m G P for all sufficiently small 6 > 0}. 

Let w be a vertex of P. One can check that the polar cone of the feasible cone of 
P at w is the cone generated by the vectors in N{A, b, v) : 

(2.1) (fcone(P,w))° = Cone(N(A, b, w)). 

Definition 2.5. Suppose P is a non-empty polyhedron defined by Ax < b. The 
recession cone of P, denoted by Kp, is 

Kp^{x\ Ax< 0}. 

2.3. Indicator functions and Brion's theorem. For a set 5* C M^. the indicator 

function [S] : M^ -5- K of S" is defined as 

We assume the readers are familiar with the definition of algebra of polyhe- 
dra/polytopes and valuation presented in [2]. The following lemma gives the two 
important equations of indicator functions of cones of polyhedra. 

Lemma 2.6 (Theorem 6.4 and Problem 6.2 in [1]). Suppose P is a non-empty 
polyhedron without lines. Then 

(2.2) [P] = y [tcone(P, w)] modulo polyhedra with lines; 

■ueVort(P) 

(2.3) [Kp] = y [fcone(P, v)] modulo polyhedra with lines, 

■ueVcrt(P) 

where Kp is the recession cone of P. 

The following lemma gives a connection between indicator functions of cones and 
indicator functions of their polars. 

Lemma 2.7 (Corollary 6.5 in [1]). Let {Ki | i e /} &e a family of cones and let Ui 
be numbers. Then 

2^0ii[Ki] = modulo polyhedra with lines 
iei 
if and only if 

2^Cii[K°] = modulo polyhedra in proper subspaces. 

Apply the above lemma to Equation (2.3) gives us the following, which will be 
useful to us later: 

Corollary 2.8. Suppose P is a non-empty polyhedron without lines. 

(2.4) [Kp] = y [(fcone(P, v))°] modulo polyhedra in proper subspaces. 

i)eVcrt(P) 
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It turns out that the multivariate generating functions define a valuation on the 
algebra of polyhedra. 

Theorem 2.9 (Theorem 3.1 and its proof in [2]). There is a map ^ which, to 
each rational polyhedron P C K'*, associates a unique rational function f{P,z) in 
D complex variables z G C^, z — (zi, . . . , zu), such that the following properties 
are satisfied: 

(i) The map ^ is a valuation. 
(ii) // P is pointed, there exists a nonempty open subset Up C C^ , such that 

X^QGPnzo 2" converges absolutely to f{P,z) for all z £ Up. 
(iii) If P is pointed, then f{P,z) satisfies 

f{p,z)= J2 ^" 

aePnz-D 

for any z £ C^ where the series converges absolutely. 
(iv) If P is not pointed, i.e., P contains a line, then f(P,z) = 0. 

Using this valuation property and Equation (2.2), we immediately have Brion's 
theorem: 

Theorem 2.10 (Brion, 1988; Lawrence, 1991). Let P be a rational polyhedron. 
Then, 

/(P,z)= Y. /(tcone(F,z;),z). 

t>eVcrt(P) 

Corollary 2.11. If P an integral polyhedron, then 

(2.5) f{P,z)= J2 zV(fconc(P,iO,z). 

■ueVGrt(P) 

Hence, for any positive integer t, 

(2.6) f{tP,z)= J2 z*V(fconc(P,«),z). 

ueVcrt(P) 

By Theorem 2.10, the problem of finding the MGF of a rational polyhedron P 
is reduced to finding the MGF of each of its tangent cone. When P is integral, the 
problem is reduced further: by Corollary 2.11, it suffices to find the formulas for 
the MGF of the feasible cone of each vertex of P. 

2.4. MGFs of unimodular cones. In general, one cannot calculate the MGF of 
a cone just by reading its generating rays. If a cone is not simple, i.e., the number 
of rays that generate the cone is larger than the dimension of the cone, one usually 
has to triangulate the cone into simple cones first. Even if a cone is simple, it is 
usually impossible to calculate its MGF directly from its generating rays. However, 
it can be done when the cone is unimodular. A pointed cone K in M^ generated by 
the rays {ri}i<i<d is unimodular if r^'s form a Z-basis of the lattice Z^ n span(/\ ). 

Lemma 2.12 (Lemma 4.1 in [2]). Suppose K is a unimodular cone generated by 
the rays {ri}i<i<d. Then 



/(A',z)=n-^ 
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Because computing the MGF of a unimodular cone is easy, it is easy to compute 
the MGF of an integral polyhedron/polytope whose feasible cones are all unimod- 
ular. Therefore, we give the following definition. 

Definition 2.13. A polytope P C M^ is totally unimodular if every vertex of P is 
a lattice point and every feasible cone of P is unimodular. 

Corollary 2.14. Suppose P C M^ is a totally unimodular polytope. Then 

«eVert(P) J=l 

where r„ i, . . . , r„^ are the generating rays of the vertex v. 

Proof. It follows from Corollary 2.11 and Lemma 2.12. D 

2.5. Obtaining volume and Ehrhart polynomial from MGF. Barvinok gives 
a polynomial algorithm to express a simple cone as signed sum of unimodular cones 
[2, Theorem 4.2]. Applying this algorithm to each simple cone in a triangulation 
of a cone, one can get the MGF of the cone. Using this idea and Theorem 2.10, 
Barvinok [2, Theorem 4.4] shows that for any P C M^ a d-dimensional rational 
polyhedron, the multivariate generating function of P is in the form of 

(2.7) /(P,z)=Ve, 



(1-Z'-'.1)---(1-Z'''.'i)' 



where ei = { — 1,1}, Vi,ri^i, . . . ,ri^d G Z^, the w^'s are all vertices (with multiple 
occurrences) of P, and Cone(ri_i, . . . ,ri^d) is unimodular, for each i. 

Because the above formula is obtained via Theorem 2.10, by Corollary 2.11, if 
P is an integral polytope, we only need to replace Vi with tvi in (2.7) to get the 
MGF of its dilation tP. Then one can use the residue calculation showed in [2] to 
find the Ehrhart polynoial and the volume of P. 

Lemma 2.15 (Lemma 5.4 in [10]). Suppose P C K^ a d-dimensional integral 
polytope and the MGF of P is in the form of (2.7). Then for any choice of c ^ R^ 
such that (c, r^ .,) ^ for each i and j, the Ehrhart polynomial of P is 

(2.8) z(P, t)^J2t^Y. H^TT r(^^' v,))hd,^ki{c, r,,i), . . . , (c, r,,,)). 



fe=0 



n,=i(c,?'jj) 



where tdd~k{ ) is a Todd polynomial. (See [2] or [10] for definition of Todd polyno- 
mials.) 

In particular, the normalized volume of P is 



(2.9) Vol(P) = 1^6 






Hence, the problem of finding formulas for the volume and Ehrhart polynomial 
of an integral polytope is reduced to finding the formula for its MGF. Therefore, 
our paper will be mainly focused on computing the MGFs of polytopes. 
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3. A PERTURBATION METHOD 

When calculating the MGF of a polytope/polyhedron which has non-simple fea- 
sible cones, we usually triangulate those non-simple feasible cones into simple cones, 
and then apply various algorithms [1, Chapter 16] for computing MGFs of simple 
cones to find the final formula. In this section, we describe a perturbation method 
that can be used to replace the triangulation step in the above procedure. The 
method is particularly useful when the perturbed polytopes are totally unimodu- 
lar, in which case instead of using other algorithms, it suffices to use Lemma 2.12 
or Corollary 2.14 to obtain the MGFs. 

Theorem 3.1. Let A be an N x D matrix and b G M.^ . Suppose P is a non-empty 
polyhedron in R^ defined by Ax < b and b(t) is a continuous function on some 
interval containing such that b(0) = b. Suppose for each t ^ in the interval, 
b(i) defines a non-empty polyhedron P{t) = {x | Ax < h{t)} with exactly £ vertices: 
Wt^i, . . . , Wt.e, and the feasible cone of P{t) at Wtj does not depend on t, that is, for 
each j '■ 1 < j < £, there exists a fixed cone Kj such that fcone(P(t), Wtj) = Kj for 
all t ^ 0. Then we have the following: 

(i) For each fixed j '■ 1 < j < i, the set of vertices {wt.j} converges to some 

vertex of P as t goes to 0. 
(ii) For each vertex v G Vert(P), let Jy be the set of j 's where {wtj} converges 

to V. Then 

[fcone(P, w)] = y [A'j] modulo polyhedra with lines. 

Therefore, 

/(fcone(P,i;),z)- Y, f{Kj,z). 

One sees that Theorem 1.1 follows directly from Theorem 3.1 and Corollary 2.11. 

The phenomenon described in Theorem 3.1 is similar to the one in Theorem 
18.3 of [1]. However, the conditions given in these two theorems are different. (For 
example, the convergence of the vertices Wt^j is a hypothesis in Theorem 18.3 of 
[1] but is a conclusion in our theorem.) One might be able to modify the proof of 
Theorem 18.3 of [1] to give a proof for part (ii) of Theorem 3.1 provided that (i) 
is proved. However, we give an independent proof in our paper. Before we prove 
Theorem 3.1, we give an example demonstrating how the theorem works and then 
give a preliminary lemma. 

Example 3.2 (Example of Theorem 3.1). Let 
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,0 <i< 1/5. 
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Ax. < b} is just the unit square with vertices 
= (0, 0), W2 - (0, 1), ^3 = (1, 1) and Vi = (1,0). 
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h = 3/20 



Figure 1 . Example of polygons converging to the unit square 



3t), 



For any t e (0, 1/5), the polygon P{t) has seven vertices 

wt,i = (t, 0), wta = (i, 1), 

Wt,3 = {l-3t,l), WtA^il-t,l-t), Wt^5^ (1,1 

wts = {l,2t), wt,7 = il~2t,0). 
See Figure 1. It shows polygons P(3/20), P(l/10) and P(l/20). One can see that 
as t decreases, the polygons look more and more similar to the unit cube P = P(0). 
In this example, as t goes to 0, 

Wt^l^Vi, Wt,2 ^ V2, 

Hence, 



W't,3,W't,4,Wt,5 -^ W3, Wt,6,Wt^7 -^ V4. 



J„, = {!}, J„, = {2}, J„3 = {3, 4, 5}, J„, = {6, 7}. 
Note that the feasible cone of P{t) at wtj does not depend on t (for t £ (0, 1/5)). 
Let Kj := fcone(P(t),u'i.j). By Theorem 3.1, we have 



(3.1) [fcone(P,ui)] 

(3.2) [fcone(F,W2)] 

(3.3) [fcone(F,V3)] 

(3.4) [fcone(P, U4)] 



[Ki], modulo polyhedra with lines; 

[K2], modulo polyhedra with lines; 

[K3] + 1X4] + [K^], modulo polyhedra with lines; 

[Kg] + [A'y], modulo polyhedra with lines. 

We remark that Equations (3.1) and (3.2) are trivial, and actually hold with direct 
equalities (without modulo equivalence). Equations (3.3) and (3.4) are less trivial. 

We need the following lemma to prove Theorem 3.1. 

Lemma 3.3. Let A be an N x D matrix with row vectors rii, . . . , njv and b G R^. 

Suppose P is a non-empty polyhedron in S.^ defined by Ax < b and h{t) is a 
continuous function on some interval containing such that b(0) = b. Let P{t) 
be the polyhedron defined by the linear inequalities Ax < h(t). Then we have the 
fallowings: 

(i) The vertex sets {Vert(P(i)) \ t =^ 0} converge to the vertex set Vert(P) of 
P; that is, for any e > 0, there exists S^, such that for any t : < |i| < Se, 
any vertex w of P{t), there exists a vertex v of P satisfying \w — v\ < e. 
(ii) Let e be strictly smaller than one half of the minimum distances between 
two vertices of P, and t : < \t\ < de, where Se is defined as in (i). For any 
vertex w of P{t), the vertex v of P satisfying \w ~ v\ < e is unique, and 
such a pair of w and v satisfies 

(3.5) N(A,b(t),u;) CN(A,b,u). 
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Hence, we have 

(3.6) (fconc(P(t),u;))° C (fconc(P,z;))°. 

(iii) Choose e and t as in (ii). Suppose P{t) is non-empty and has vertices. For 
any vertex v G Vcrt(P), let 

Wt,v = {we YcTt{P{t)) I \w - v\ < e} 

be the set of vertices of P{t) that are within e-distance to v. Then 

[iconc{P,v)] = y [iconc{P{t),w)] modulo polyhedra with lines. 

Therefore, 

/(fcone(P,«),z)= Y. /(fcone(P(i),«^),z). 

Proof. Any D linearly independent rows of A form an invertible matrix. Let 
Ml, . . . , Mk be all the invertible matrices arising this way. For each j : 1 < j < k, 
let 

Ij ~ { the indices of the rows from A that form matrix A/j}, 
and let b-'(i) = {bi{t))i^j. be the corresponding sub- vector of h{t) and let b^ = 
{bi)i(=i- be the corresponding sub-vector of b. Note that b^(i) is a continuous func- 
tion where b-' (0) = b-' . 

Let e > 0. For any j : I < j < k, because AP"^(b^(i)) is a continuous function 
on t, there exists Sj such that 

(3.7) 0<\t\<6, => \Mr\hn-M-\h^m<e. 

If M^ (b-*) is not a vertex of P, there exists i ^ Ij such that ni{M^ (b-')) > bi. 
Because ni{M~^(h^ (t))) is a continuous function on t, there exists (5' such that 

(3.8) < |t| < 5'^ => n,(A/ri(b^(t))) > b,. 
We choose 

4 = mm{{5, I 1 < .? < k}, {S^ I A/ri(b^) ^ Vert(P)}). 

We claim for any t : < |i| < S^, and any vertex w of P{t), there exists a vertex 
u of P such that |u) — w| < e. After all, since w is a vertex of P{t), for some (not 
necessarily unique) j, we have that 

(3.9) w = Ar\h^t)) 



(3.10) n,(A/ri(bJ)) < 6„ Vi ^ /,. 



Hence, A'/-"^(b^) has to be a vertex of P because otherwise J' exists and (3.8) 
contradicts with (3.10). Then our claim, i.e., (i), immediately follows from (3.7). 

We now prove (ii). The uniqueness choice off follows from the triangle inequality. 
Suppose ni(x) = bi{t) is a supporting hyperplane for w. Then 3 j : 1 < j < fc such 
that Ui is a row vector oi Mj and w ~ AI^^{h^{t)). By our argument above, M~^{h) 
is the unique vertex v oi P that is within e-distance to w. Hence, ni(x) = 6^ is a 
supporting hyperplane for v = Af~^(b). Therefore, (3.5) follows. The inclusion 
relation (3.6) follows from (3.5) and (2.1). 
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Finally, we prove (iii). Suppose P{t) is non-empty and has vertices. It is clear 
that P is a polyhedron with vertices, thus without lines. Also, P and P{t) have 
the same recession cone: 

K ■.= {xeR^ \ Ax< 0}. 

By (3.6), we have 

(fcone(P(t),u;))° C (fconc(P, w))°, Vw e Wt,v 

Thus, 

IJ {iconc{P{t),w)y C (fcone(P,i;))°. 

However, applying (2.4) to both P and P{t), we get 

J2 [(fcone(P,«))°]^ [A-] ^ E E [(fcone(P(i),u;)1 

■ueVert(_P) veYeTt{P) wGWt,^ 

modulo polyhedra in proper subspaces. 
Hence, we must have 
[(fcone(P, w))°] = 2_. [(fcone(P(i), u;)°] modulo polyhedra in proper subspaces. 

w£Wt,v 

Then (iii) follows from Lemma 2.7, the above formula, and Theorem 2.9. D 

Proof of Theorem 3.1. Clearly, (ii) follows from (i) and Lemma 3.3/(iii). So we just 
need to prove (i). 

Fix j G {!,... ,^}. By Lemma 3.3, for sufficiently small e, there exists 5^, such 
that for any t : < |i| < 8^, any Wt^j is within e-distance to a unique vertex, say 
Wt, of P. It suffices to show that vt are the same vertex for any t : < |t| < Jg. We 
assume to the contrary that there exist ti ^ ^2 such that vt-^ ^ vt^ . However, by 
(3.6), we have 

K° C (fcone(P, wt J)° and K° C (fcone(P, WtJ)°. 

Hence, 

K° C (fcone(P,wtJ)° n (fcone(P, wtj)°. 
Since Kj is a pointed cone, the polar K° is full-dimensional. However, the intersec- 
tion of (fcone(P, Vt^))° and (fcone(P, utj))" lies in the hyperplane {n | n{vt^ ~"t2) = 
0}. This is a contradiction. Hence, we have (i). D 

Remark 3.4. Theorem 3.1 and Theorem 1.1 still hold if we replace the continuous 
function b(t) with a sequence {bi, b2, . . . } converging to b. 

4. Properties of Transportation Polytopes 

In Section 5, we will apply the perturbation method introduced in the last section 
to transportation polytopes. Before that, we review properties of transportation 
polytopes. Most of them can be found in [13]. 

Transportation polytopes have natural connection to the complete bipartite 
graphs. Let Km.n be the complete bipartite graph with m vertices ui, . . . ,Um 
on the left and n vertices wi, . . . ,Wn on the right. Throughout the rest of the 
paper, we refer to Mi's as the left vertices and Wj's as the right vertices. Denote by 
Cij the edge in Km,n connecting Ui and the Wj. For any subgraph G of K,n,m we 
denote by E{G) the edge set of G. 



12 FU LIU 

Let Am,n be the {m + n) x mn incidence matrix of Km.n- (Then the column 
AJ;f „ of A„i^n corresponding to the edge e^j is the (?7i + n)-vector where the ith and 
(?7i + j)th component arc 1 and zero elsewhere.) Then the transportation polytope 
T(r, c) can be described by 

A-m,n^ — T \ ^ X > 0. 

yc / 

We often call matrix A^.n the constraint matrix of the transportation polytopcs of 
order m x n. 

A vertex of a transportation polytope of order m x ri is non-degenerate if it 
has exactly m + n ~ 1 entries that arc positive; otherwise it is degenerate. A 
transportation polytope is non- degenerate if all its vertices are non-degenerate; 
otherwise it is degenerate. It's easy to see that every non-degenerate transportation 
polytope is simple. It is known that Am,n is a totally unimodular matrix^ i.e., every 
minor of A,„_„ is 0, 1, or —1. Therefore, if a transportation polytope is non- 
degenerate or simple, it is a totally unimodular polytope. Then we can apply 
CoroUary 2.14 to find its MGF. 

There is an easy way to identify non-degencratc transportation polytopcs. 

Theorem 4.1 (Theorem 1.2 of Chapter 6 in [13]). The transportation polytope 
T(r,c) is non- degenerate if and only if the only nonempty index subsets I C [m] 
and J C [n] satisfying X^ie/ ''« = Sfej '^j '^^^ ^ ~ ["^] '^"-'^ J ~ ["]• 

To express the MGF of a totally unimodular polytope, wc need to know how to 
describe its vertices, and the generating rays of each feasible cone. 

4.1. Vertices of transportation polytopes. The vertices of a transportation 
polytope can be characterized with the auxiliary graphs. 

Definition 4.2. The auxiliary graph of M, denoted by aux(M), is the induced 
subgraph of A'„i_„ with edge set {eij \ M{i,j) > 0}. 

We denote by vertAux(r, c) the set of all the auxiliary graphs obtained from 
vertices of T{r, c). 

A point M G 7~(r, c) is a vertex of T(r, c) if and only if the column vectors 
{^m,n I M{hj) > 0} are linearly independent. However, using matroid theory, one 
can show that given any index set ind C [m] x [n] , the column vectors {Al^^ \{hj) G 
ind} are linearly independent if and only if the induced subgraph of K^^ with edge 
set {cij- I {i,j) S ind} is a forest. Therefore, we have the following theorem and 
lemma, both of which follow from Theorem 2.2 of Chapter 6 in [13]. 

Theorem 4.3. Let M G T(r, c). Then M is a vertex of T(r, c) if and only if 
aux(M) is a spanning forest of K„i,n- Furthermore, aux induces a bijection between 
the vertex set of T(r, c) and the set of spanning forests of Km.n that are auxiliary 
graphs of some points in T(r, c). 

In particular, if T{r, c) is non- degenerate, M is a vertex of 7~(r, c) if and only 
if aux(M) is a spanning tree of Km,n- Thus, aux induces a bijection between the 
vertex set o/T(r, c) and the set of spanning trees of Km,n that are auxiliary graphs 
of some points in T{v, c). 

Lemma 4.4. Let M and N be two vertices in 'T{v, c). Suppose there is a spanning 
tree T of Km,n such that both aux(M) and aux(A^) are subgraphs of T. Then M 
and N are the same vertices. 
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We remark that a weak statement of Theorem 4.3 is that aux induces a bijec- 
tion between Vert(7~(r, c)) and vertAux(r,c). This is the first reason why the set 
vertAux(r, c) is important (and why we have defined this set). In fact, in addition 
to the one-to-one correspondence to the vertex set of T(r, c), the set vertAux(r, c) 
also determines the combinatorial structure of T(r, c) as we will see in the next 
subsection (cf. Corollary 4.9). 

4.2. Feasible cones of transportation polytopes. Section 4.1 of Chapter 6 in 
[13] gives a complete description of the generating rays of feasible cones of trans- 
portation polytopes, as well as a characterization of when two vertices are adjacent 
in a transportation polytope. We summarize the results as two lemmas below 
(Lemma 4.6 and Lemma 4.7). We begin with a preliminary definition. 

Definition 4.5. Let T be a spanning forest of K„in and e ^ E{T). Suppose T U e 
creates a cycle. Note this cycle is unique if it exists. The edge e must be contained 
in the cycle. Suppose this cycle is: 

We define cycle(r, e) to be the m x n matrix whose entries are defined as 

1, if («,.?) e {{ii,ji), {i2,.h), ■ • • , (Js,js)}; 

(4.1) cycle(r, e)(z, j) = <( -1, if (*, j) e {{t^Ji), (13,^2), ■ • ■ , (ii, Js)}; 

0, otherwise. 

Lemma 4.6. Let M be a vertex in 'T{r,c), and T = aux(M) the auxiliary graph 
of M. (By Theorem 4-3, T is a spanning forest of Kmn-) Then 

(4.2) {cycle(r, e) | T U e creates a cycle} 

is the set of rays that generates the feasible cone of T{v^ c) at the vertex M. 

In particular, i/7~(r,c) is non-degenerate. Then the generating set (4.2) can be 
rewritten as 

(4.3) {cyclc(T,e) I e^£;(T)}. 

Lemma 4.7. Let M and N be two distinct vertices in 7~(r, c). Then Ad and N are 
adjacent if and only if the union of aux(M) and aux(7V) has a unique cycle. 

Moreover, if M and N are two adjacent vertices, the unique cycle in the union 
o/aux(M) and aux(iV) is the same as cycle(aux(M), e), for some e ii^(aux(M)). 
Thus, the unique cycle determines the ray from M to N as described in (4.1). 

It turns out that the transportation polytope T(r, c) is integral if r and c are both 
integer vectors. Hence, the following corollary follows immediately from Lemma 
2.12, Corollary 2.14 and Lemma 4.6. 

Corollary 4.8. Suppose T(r,c) is a non-degenerate transportation polytope. For 
any vertex M o/7"(r, c), we have 

(4.4) /(fcone(r(r,c),Af),z)= J] 1 _ ,cyci!(aux(M).e) ■ 

e^B(aux(j\/)) 

Suppose further that r and c are integer vectors. Then the multivariate generat- 
ing function of T(r, c) is 

(4.5) /(r(r,c),z)= Y. ^"' n 1 _ ,cycleiux(M).e) ■ 

A/eVcrt(r(r,c)) e0_E(aux(A/)) 
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We also have the fohowing corollary to Lemma 4.7. 

Corollary 4.9. Suppose 7~(r,c) and 7~(r',c') are two transportation polytopes of 
order m x n satisfying 

vertAux(r, c) = vertAux(r',c'). 

Let T e vertAux(r, c), and Mt and M^ be the vertices o/ 7~(r, c) and T{r',c'), 
respectively, corresponding to T . Then 

fcone(r(r, c), A/t) = fcone(r(r', c'), A^). 

5. A PERTURBATION FOR TRANSPORTATION POLYTOPES 

We describe a perturbation that works for any transportation polytope. 

Lemma 5.1. Suppose r = (ri, r2, . . . , rm) and c — (ci, C2, . . . , c„) are two rational 
positive vectors. We define 

1 



r(i) = (ri -i, ...,r„i -i), c(t) = (ci, . . . ,c„_i, c„ - mi), 0<i< 



Km 

where K is the greatest common divisor of the denominators of r and c. Then we 
have the following: 

(i) For any t G (0, tt—), the transportation polytope T{v{t), c(i)) is non- degenerate. 
(ii) For any t £ (0, jy — ), the set vcrtAux(r(t), c(t)) is independent oft. 
(iii) {T(r(i),c(i)) | t G (Oi"^^)} ^^ a family of transportation polytopes satisfy- 
ing the condition of Theorem 3.1. 

Remark 5.2. The perturbation we defined in Lemma 5.1 perturbs the original poly- 
tope by i(-l, -1, . . . , -1, 0, . . . , 0, -m)^. The vector (-1, -1, . . . , -1, 0, . . . , 0, -m)^ 
is in the cone spanned by the column vectors of the constraint matrix Am.n- We 
expect that for any generic vector p in the cone spanned by the columns of Am.n, 
perturbing the original polytope by tp works. However, we choose this particular 
one since it gives nice combinatorial results for the central transportation polytopes 
of order kn x n, which we will discuss in the next section. 

We give the following two definitions before the proof of the above lemma. 

Definition 5.3. Let a; be a real number. The floor or integer part of x, denoted 
by [xj is the biggest integer that is not greater than a;, and the fractional part of x, 
denoted by frac(x), is x— \_x\ . The ceiling of x, denoted by [x] is the smallest integer 
that is not smaller than x, and the co-fractional part of x, denoted by cofrac(x), is 

[x] — X. 

Definition 5.4. Let G be a subgraph of the complete bipartite graph K„i^n- 
Let dj be the number of edges in G connecting to Wj , the jth right vertex of Km^n- 

We caU (di, (i2, • ■ • , dn) the right degree sequence of G and write S{G) = {di, . . . , dn). 
We also define 

/(G) :~ the number of u^'s in G, i.e., the number of left vertices in G. 

For convenience, for any spanning tree T of Km,n, we consider T as a rooted 
tree rooted at Wn, the nth right vertex. For any vertex v of T, we denote by T„ the 
subtree of T rooted at v. 
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Proof of Lemma 5.1. Because T{Kr{t),Kc{t)) is a dilation of 7~(r(t), c(i)), the 
polytopes T{Kr{t),Kc{t)) and T{r{t),c{t)) have exactly the same combinatorial 
structures. Therefore, without loss of generality, we can assume that r and c are 
integer vectors, and K = 1. 

(i) Suppose 7^ / C [m] and 7^ J C [n] are two index sets satisfying 

(5.1) 5^r(i),=^c(t),. 

lei je.J 

The co-fractional part of the left hand side of (5.1) is \I\t ^ 0. Hence, n ^ J, 
because otherwise the right hand side of (5.1) is an integer. Then the co- 
fractional part of the right hand side of (5.1) is mt. Therefore, |/| = m 
and / = [m]. This implies that J = [n]. Therefore, by Theorem 4.1, the 
polytope T(r(i),c(i)) is non-degenerate, 
(ii) Let to e (0, -^ = ^), and let T e vertAux(r(io),c(to))- By Theorem 4.3, 
T is a spanning tree of Km^n- It suffices to show that for any t e (0, — ), 
the tree T is the auxiliary graph of a vertex of T(r(i), c(t)). Let Mxito) be 
the vertex of 7~(r(to),c(to)) corresponding to the tree T. 
We claim that, considering T a rooted tree rooted at «;„, 

cofrac(AfT(^o)(*5 j)) ~ K'^ui) ^Oi if Wj is the parent of Ui in T; 
frac(Af7^(to)(*, j)) ~ l{Tw-) to, if Ui is the parent of Wj in T; 

and the matrix Mt whose entries defined by the following equation is a 
vertex of T(r, c). 

{\MT{to){i,j)~\, if Wj is the parent of Ui in T; 
[MT{to){'i,j)\, if Ui is the parent of Wj in T; 
0, otherwise. 

The claim can be proved by induction on hook lengths of vertices of T. 
(Recall the hook length of a vertex i; in a rooted tree is the number of 
descendants of v.) We define the matrix MT{t) as follows: 

{^T{hj) — IC^Ui) t, if Wj is the parent of Ui in T; 
MriiJ) + l{Twj) t, if u^ is the parent of Wj in T; 
0, otherwise. 

It is clear that Mxit) is a vertex of T(r(i),c(i)) with auxiliary graph T. 

See Figure 2 for an example. The tree on the right side of the figure is 
the same one as the one on the left; it is just drawn as it is a rooted tree 
rooted at ^3. The solid lines are the edges where some Wj is the parent and 
the dashed lines are the edges where some Ui is the parent. The number 
next to an edge {u,;, Wj} is the entry of MT{t){i,j). 
(iii) It follows from (ii) and Corollary 4.9. 

D 

It turns out we can use auxiliary graphs to determine which vertices of the 
perturbed polytopes T(r(i),c(i)) (described in Lemma 5.1) converge to a given 
vertex of T(r,c). 
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(ri + • • ■ + rs - ci - C2)-/U\'6 - * 

! (ci + C2 - ri - ra - rg - rg) + 4t 

= n - t / \(r2 + r-a + rg - c2) - 3t 

1^2 



Ul 




I (c2-r3-r5) + 2i 

Figure 2 

Lemma 5.5. Assume the conditions of Lemma 5.1. Let t^ 6 (0, -J— ), a^ii tet 
T G vertAux(r(to), c(to))- Define Mt and Mxit) as in (5.2) and (5.3), respectively. 
Then 

(i) linit^o MT(t) = MT. 

(ii) Mt is the unique vertex of 7"(r, c) satisfying that its auxiliary graph is a 
subgraph of T. 

Proof. It is clear that Mxit) converges to Mt, and Mt is a vertex of 7~(r, c) 
satisfying that its auxihary graph is contained in T. It is left to show the uniqueness. 
However, the uniqueness follows from Lemma 4.4. D 

Lemma 5.5 tells us which vertices of the perturbed transportation polytopes 
T(r(i),c(i)) converge to a given vertex Af of the original transportation polytope 
r(r,c): 

lim Afy(i) = M <=^ aux(M) is a subgraph of T. 

Now we can describe the MGF of the feasible cone of each vertex of T(r, c) 

Corollary 5.6. Assume the conditions of Lemma 5.L Let t e (0,-^— ). For any 
vertex M o/T(r,c), let 

(5.4) PertAux(M) := {T G vertAux(r(t), c(i)) | aux(A/) is a subgraph ofT} 

be the set of auxiliary graphs of vertices of the perturbed transportation polytope 
T{r(t),c{t)) that contain the auxiliary graph of AI as a subgraph. 
Then 

(5.5) [fcone(r(r,c),M)] ee ^ [fcone(r(r(i), c(i)), AfT(i))] 

T6PertAux(A/) 

modulo polyhedra with lines, 
where MT{t) is defined as in (5.3). 
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Hence, 

(5.6) /(fcone(r(r,c),Af),z)= ^ \{ ^ ^ ^^,,,^r,e) - 

TePcrtAux(J\f ) e<^E{T) 

Proof. Formula (5.5) follows from Lemma 5.1/(iii), Theorem 3.1, and Lemma 5.5. 
Formula (5.6) follows from (5.5), Lemma 5.1/(1), and Formula (4.4). D 

Wc can also describe the MGF of an integral transportation polytope. 

Corollary 5.7. Assume the conditions of Lemma 5.1 and further assume that r 
and c are integer vectors. Let t G (0, — ). Then the multivariate generating function 
ofT{r,c) is 

(5.7) /(r(r,c),z)= Y: -"^ n 1 _ ,c'cIe(T.e) ^ 

TevortAux(r(t),c(t)) e^E{T) 

where Mt is defined as in (5.2), or equivalently, Mt is the unique vertex ofT{v,c) 
satisfying that its auxiliary graph is a subgraph of T. 

Proof. Formula (5.7) follows immediately from Lemma 5.1, Theorem 1.1, and Lemma 
5.5. D 

Maximum number of vertices. We finish this section by an additional result 
we obtain from the perturbation we define in Lemma 5.1. 

Suppose T{y, c) is a central transportation polytope of order mxn, and T(r', c') 
is a transportation polytope of same order. By Theorem 7.1 of Chapter 6 in [13], 
T(r', c') has the maximum possible number of vertices if and only if for all A € (0, 1), 
the transportation polytope 7"(Ar + (1 — A)r', Ac + (1 — A)c') is non-degenerate. 
Because the perturbation we define in Lemma 5.1 is linear, we have the following 
result: 

Lemma 5.8. Suppose 7~(r, c) is a central transportation polytope. Then the trans- 
portation polytopes T{r{t), c{t)), < t < -j^, we defined in Lemma 5.1 achieve the 
maximum number of vertices among all the transportation polytopes of order mxn. 

Proof. For any A G (0, 1), it is clear that 

r(Ar + (1 - A)r(t), Ac + (1 - A)c(t)) = r(r((l - A)t), c((l - \)t)). 

Since < (1 — \)t < t < -j^, by Lemma 5.1/(i), the above transportation polytope 
is non-degenerate. Then the result follows. D 

6. Central transportation polytope of order kn x n 

In this section. We always assume that 7~(r, c) is a central transportation poly- 
tope of order mxn, where m = kn, and r = (a, . . . , a) and c = (6, . . . , &) are two 
integer vectors. We apply the perturbation defined in the last section to family of 
central transportation polytopes {T(r(i), c(i)) \ t G (0, — )} satisfying the condition 
of Theorem 3.1, where 

r(t) = (a - t, . . . , a - t), c(t) ^ {b,...,b,b- mt), < t < — . 

771 

We have the following theorem on the vertices of T{r{t),c{t)). Recall that STk,n 
is the set of spanning trees of Kkn,n satisfying 5{T) = {k + 1, . . . ,k + \,k). 
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Theorem 6.1. Let t G (0,—). The set of vertices of T{r{t),c{t)) is in bijection 
with the setSTk,n- More specifically, STk.n is the set of the auxiliary graphs of the 
vertices o/7~(r(t), c(t) : 

STk,n = vertAux(r(i),c(t)). 

Furthermore, for any T S STk,m the corresponding vertex of 7'{r,c) is the matrix 
Mxit) whose entries are defined as below, considering T a tree rooted at Wn-' 

{a — l{Tui) t, if Wj is the parent of Ui in T; 
l{Tw-) t, if Ui is the parent of Wj in T; 

0, otherwise. 

Proof. Wc have akn — bn. So ak — b. 

Suppose M is a vertex of T{r{t),c{t)). By Theorem 4.3, the auxihary graph 
T := aux(Af ) is a spanning tree. Hence, the number of edges in T is kn + n ~ 
1 = (fc + l)n — 1. It is clear that each entry of M cannot exceed a ~ t. For any 

J : 1 < .7 < ^ ^ li since X)i^i-^(*'i) ^ ^ ^^^'^ i^ ^ 0^ < ^i there arc at least 
k + 1 Vs such that M{i,j) is positive. This means the right vertex Wj has at least 
k + 1 adjacent edges in T. Similarly, we can argue that the right vertex w„ has 
at least k adjacent edges. However, the sum of the right degree sequence of T is 
equal to the number of edges in T which is (fc + l)n — 1. Thus, we must have 
S{T) = {k + l,...,k + l,k). 

On the other hand, given T G STk,m to show that T is the auxiliary graph of a 
vertex of T(r, c), it suffices to verify Mt defined by (6.1) is a point in 7"(r, c). This 
can be proved directly by checking the row sum of Mt is always a and the column 
sum is always b = ak. D 

We see that Theorem 1.2 follows from Theorem 6.1 and Corollaries 5.6 and 5.7. 
We can actually analyze the data used to enumerate the formulas in Theorem 1.2 
further, which helps to restate Theorem 1.2 with more fundamental combinatorial 
objects, as well as to figure out the number of vertices in T{r{t),c{t)). It is clear 
that the vertices of T(r, c) are the {0, a}-matrices in which each row has exact one 
entry of a and each column has exactly k entries of a. These matrices corresponding 
to "fc to 1" matching from the m = kn left vertices to the n right vertices. This 
motivates the following definition. 

Definition 6.2. We call an kn x n matrix M a k-to-1 matching matrix if M is a 
{0, l}-matrix such that there is exactly one 1 in each row and exactly k I's in each 
column. We denote by Aiatk,n the set of all the kn x n fc-to-1 matching matrices. 
We also call the auxiliary graph of each fc-to-1 matching matrix a k-to-1 matching 
graph. 

With this definition, the vertices of T(r, c) is the set 

Vert(r(r,c)) == {aM \ M £ Matk,n} =: aMatk.n- 

Now we connect the set vertAux(r, c) = STk,n with simpler combinatorial ob- 
jects, denoting by TZn the set of all the rooted trees on {'Wi,W2, . . . , w„} rooted at 

Wn. 

Lemma 6.3. There is a bijection between AAatk^n x Ti-n x [fc]"~^ and STk,n- 
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Proof. Given M e Matk.n, a rooted tree R e Tin and f = (/i, . . . , /„-i) £ [A:]""^, 
we can eonstruct a tree T E STk,n from (Af , i?, f ) in the following way: We start 
with the fc-to-1 matching graph aux(M). If Wj^ is the parent of Wj in _R, we add an 
edge connecting Wj and the (/j)th left vertex that is matched to Wj^ in aux(A'/). 
After adding these n — 1 edges, one can check that we actually obtain a spanning 
tree T of Kkn.n rooted at w„ and each right vertex Wj has exactly k children. We 
can ignore the root. Then T is in STk.n- 

For example, the tree in Figure 2 is the image of the tuple (M, R, f ) shown in 
Figure 3. 



M 



(I 0\ 






1 


• ^3 


1 
1 


-^ = OWi 


1 
\q l) 


• W2 




Figure 3 





f =fl,2 



Conversely, let T G ST^^- Considering T a rooted tree rooted at w„, because 
5{T) = (fc + 1, . . . , fc + 1, fc), one sees that each right vertex Wj has exactly k children. 
We separate the edge sets of T into two sets: 



El 

E2 



{e e E{T) 
{e G E{T) 



e connects a right vertex Wj to one of its children}; 
e connects a right vertex Wj to its parent}. 



In the tree shown in Figure 2, Ei is the set of solid lines and E2 is the set of the 
dashed lines. 

It is clear that Ei is the edge sets of a k-to-l matching graph. Let M be the 
fc-to-1 matching matrix corresponding to it. If we contract all the edges in Ei in T 
(and remove all the left vertices Ui), we obtain a rooted tree R G TZn- Finally, we 
define f G [/c]""-^ from E2: If {wj,Ui} G E2, i.e., Ui is the parent of Wj in T, and 
suppose Ui is the sth child for its parent, we define fj = s. Hence, we defined the 
inverse map. D 

We denote by <& the bijection from Matk^n x Tin x [fc]"~^ to STk^n given in the 
above proof. Let w be a vertex of 7~(r, c), i.e., v = aM for some M G Aiat^.n- 
Recall that the set PertAux(t;) given in Theorem 1.2 is defined as 

PertAux(w) = {T G STk,n \ the auxihary graph of w is a subgraph of T} 
= {r G STk,n I the auxiliary graph of M is a subgraph of T.} 

It is clear that for any fixed vertex v = aM of T(r, c), the bijection $ induces a 
bijection from 7^„ x [fc]"-i to PertAux(i;), which maps (T, f ) to $(A/,r,f). We 
denote this map by $a/. Therefore, 

PertAux(f) = $A/(7^„ x [fc]""^). 

We now restate Theorem 1.2. 
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Corollary 6.4. Suppose M g Matk,n- Then the MGF of the feasible cone of 
T(r, c) at V = aM is 

(6.2) /(fcone(r(r,c),z;),z)= Y. IT i _ ^cycie(T,e) ■ 
We have two formulas for the MGF o/T(r,c); 

(6.3) /(r(r,c),z)= Y. -""' E n 1 „ ,c'cle(T,e) ^ 

and 

(6.4) /(r(r,c),z)^ Y -"'^ n 1 _ ,c'cle(T,e) ^ 

TGSTfe,^ e^_E(T) 

where Mt is the matrix whose corresponding graph is the unique k-to-l matching 
subgraph ofT, or equivalently, Mt is the k-to-l matching matrix entry in the tuple 
$-i(T). 

Proof. Formulas (6.2) and (6.4) follows from Theorem 1.2 and the discussion above. 
Formula (6.3) follows from Formula (6.2) and Corollary 2.11. D 

When k ~ I and a = 1, the set A4k.n is actually the symmetric group &n, and the 
polytope T(r, c) is the Birkhoff polytope i?„. Therefore, for any a € &n, the map 
<I>cr is a bijcction from TZn to PertAux((T). Hence, we obtain the following theorem 
for the Birkhoff polytopes, which is equivalent to Theorem 1.1 and Corollary 4.1 in 

[10]. 

Theorem 6.5. Suppose a £ S„ is a vertex of the Birkhoff polytope Bn. Then the 
MGF of the feasible cone of Bn at a is 



/(fcone(B„,a),z)= ^ [| ^— ^ 



cycle(T,e) ' 

Te$„(R„)e^_B(T) 



We have two formulas for the MGF of Bn '■ 



cycle(T,e) ' 
<TGe„ Te*„(R„) e^£;(T) ~ 



/(B„,z). Y ^"^ n 13^ 

TG5Ti.„ eiE(T) 

where Mt is the matrix whose corresponding graph is the unique \-to-\ matching 
subgraph of T. 

Next, we give formulas for the number of vertices of T(r(i), c(<)). 

{kn)\ 



Corollary 6.6. The number of vertices o/ 7~(r(0), r(0)) — T(r,c) is 



{kiy 



(knV 
The number of veHices o/r(r(t), c(t)) (t e (0, ^); is TrjNTr""'^'^""^- 
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Proof. The first statement follows from the observation that the cardinality of 

It is well-known that the number of rooted trees on n vertices with a fixed root 
is n"^^. Then the second statement follows from Lemma 6.3. D 

By Lemma 5.8, we obtain another known result. 

Corollary 6.7 (Corollary 8.6 of Chapter 6 in [13]). The maximum num,ber of 

(kn)\ 
vertices among all the transportation polytopes of order kn x n is n"~ fc"~ . 

Finally, we apply Lemma 2.15 to Formula (6.4) to obtain formulas for the volume 
and Ehrhart polynomial of the central transportation polytope of order kn x n. 

Corollary 6.8. Let C be a knxn matrix such that (C, cycle(T, e)) ^ for any pair 
ofTE STk^n o,nd e ^ E{T). Then the coefficient ofV in the Ehrhart polynomial of 
the central transportation polytope 7~(r, c) of order kn x n is given by 

1 ^p ((C,aA/r))'td(fc„_i)(„_i),,((ficycle(r,e)) | e ^ EjT)) 
*^J^^„ Ue^EiT){C,cyde{T,e)} 

In particular, the normalized volume of T(r, c) is given by 

v°i(r(r,c)) = — -— -— Y. 



{{kn - l)(n - 1))! ^^ Yle^EiT) {C, cycle(r, e)) ' 



We remark that when fc = 1, the above corollary gives an equivalent result to 
Corollary 1.2 in [10]. 
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